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Abstract 

We discover quantum Hall like jumps in the saturation spectral rigidity in the semiclassical 
spectrum of a modified Kepler problem as a function of the interval center. These jumps correspond 
to integer decreases of the radial winding numbers in classical periodic motion. We also discover 
and explain single harmonic dominated oscillations of the level number variance with the width 
of the energy interval. The level number variance becomes effectively zero for the interval widths 
defined by the frequency of the shortest periodic orbit. This signifies that there are virtually no 
variations from sample to sample in the number of levels on such intervals. 
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I. INTRODUCTION 



Level correlations in the semiclassical spectra of classically integrable systems have re- 
cently received a renewed attention. The most important development was the realization of 
the long-range nature of such correlations which was explored for rectangular billiards. 
While it had been previously known that the short-range correlations are absent, the fact 
reflected by the Poisson statistics of the nearest neighbor level spacings the evidence 
for the long-range correlations was indirect, namely, through the saturation property of the 
spectral rigidity In [l| the correlation function of the level density was obtained, 

which explicitly describes the long-range correlations in the energy spectrum. Furthermore, 
in terms of an easily measured quantity, the variance of the number of levels on an energy 
interval was investigated and was shown to have very unusual properties. Namely, for the 
interval width narrower than the energy scale associated with the inverse time of the short- 
est periodic orbit (traversal along the smaller side of the rectangle), the variance equals, 
in the lowest approximation, to the mean number of levels in the interval, indicating the 
absence of correlations in level positions. For intervals wider than such energy scale, the 
variance exhibits non-decaying oscillations around the "saturation value" with the ampli- 
tude smaller, yet parametrically comparable to the latter and with the "period" of the same 
order as the above mentioned scale (the width of the interval at which the transition from 
the uncorrelated to correlated behavior occurs). 

While such behavior of the variance had being previously predicted via a formal math- 

n n 

ematical approach [5j, ref. pj established that it is a direct consequence of the long-range 
correlations between energy levels. Two independent analytical derivations were produced 
jj: one based on the direct use of quantum mechanical expressions for the energy levels 
for a particle in a box and the other based on semiclassical periodic orbit theory {3]. 
Within the latter, it was shown that the oscillations of the variance can be explained by 
just a few shortest periodic orbits. These results were confirmed numerically with the use 
of an ensemble averaging procedure wherein the rectangles of the same area, but varying 
aspect ratios, were used. The reason why the oscillations of the variance can be considered 
counter-intuitive is because with the increase of the interval width, and the corresponding 
increase of the mean number of levels, the fluctuation of the number of levels in the interval 
may actually decrease. 
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It should be pointed out that the only reason that the variance does not become zero 
for a rectangular box is that the harmonics that correspond to the shortest periodic orbits 
have incommensurate frequencies and thus add incoherently |]| . If one could find a system 
where variance would be dominated by a single periodic orbit harmonic, it could be near 
zero for certain widths of the energy interval. This would be even more counter-intuitive 
as statistically independent systems would produce different level structures yet the total 
number of levels for such intervals would be nearly the same! In this work we report that we 
found just such a system - modified Coulomb problem - and, in addition to examining the 
variance, we also show that the saturation value of the spectral rigidity [7] exhibit quantum 
jumps associated with the change in the winding number ratio of radial and angular motions 
of periodic orbits. 

In what follows, we first discuss level correlations in the semiclassical spectrum and illus- 
trate it by a classically integrable problem of a particle in a rectangular box Next, we 
give a detailed analytical description of the modified Coulomb problem and derive expres- 
sions for the saturation spectral rigidity and for the level number variance using the periodic 
orbit theory. We then proceed with the numerical evaluation of these quantities for a model 
spectrum, that captures key features of our model, where we observe the quantum jumps in 
saturation rigidity and single-harmonic oscillations of the variance. 



II. PERIODIC ORBIT THEORY OF LEVEL CORRELATIONS 

We will consider intervals [e — E/2, e + E/2] , E <C e, where the states with energies 
near e have large quantum numbers and can be described semiclassically. Denote by M (e) 
the cumulative number of levels (or spectral staircase) 3] 

Af(s) = J2Hs-e k ) (1) 

k 

where 9 is unit step function and k labels the energy eigenstates. A 'universal' (flattened) 
representation of the ladder data is obtained by rescaling the energy variable so as to elim- 
inate the particular shape of the average (J\f (s)) from the ladder^. To do this, define the 
new scaled dimensionless energy variable e' by 

s^e'{e) = (N{e)) (2) 
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Here () denotes the ensemble average. In particular, to a computed eigenvalue Sk the value 
e' k = (jV(£fc)) is assigned. Since e' is a monotone function of e 



(Af'(e>)) = (Af (e))=e> 



(3) 



so that the mean level density (and the mean level spacing - with caveats explained in Ref. 
is unity in the scaled variable 



<p'(^)> = (E^' 

\ k 

a = (py 1 = 1 



(4) 
(5) 



Since (Af (e)) is used as the scaled energy axis variable in order to give the ladders an ap- 
proximately 45° slope, it is important to have a fair idea of its functional form for presenting 
the numerical data. In what follows, we will omit the "prime" for the scaled energy variable. 

Since the present work concentrates on the long-range correlations in the spectrum of the 
eigenvalues, we will not discuss per se the distribution of nearest-neighbor level spacings, 
except to state that our computations for the modified Coulomb problem below give the 
Poisson statistics Q], as ought to be the case for a classically integrable system. For the 
statistics of large numbers of levels, the following standard measures will be used. The first 
is the spectral rigidity A 3 , defined in as the best linear fit to the spectral staircase in 

the interval [e-E/2,e + E/2] 

min 1 



A 3 (e;E)-- 
the explicit form of which is 

I r e+E/2 



(A,B)Ej e „ E/2 



de[Af(e) - A -Be}' 



(6) 



E 



deAf 2 



e-E/2 



1 

£2 



e+E/2 



deAf I 



e-E/2 



12 

£4 



e+E/2 



-1 2\ 



deeAf (e) 



e-E/2 



(7) 



For the number of levels iV on the interval [e — E/2, e + E/2] 



N(e;E)=Af(e+^) 



Afle-§ 



(8) 



the variance 

Z(e;E) = ((N-(N)f) (9) 

is another measure of the fluctuations. Notice that in the flattened spectrum ([5]), considered 
here, (N) = E. 



The fluctuation measures S and A3 can be expressed in terms of the correlation function 
of the density of levels, 

K(e 1 ,e 2 ) = (5p(e 1 )5p(e 2 )} (10) 
5p(e)=p(e)-(p(e)) (11) 

regardless of the form of K (ei, e 2 ) , for instance, 

re+E/2 re+E/2 

E{e;E)= / K{e 1 ,e 2 )de 1 de 2 (12) 

Je-E/2 Je-E/2 

Using these relationships one can further show that E supersedes A 3 via an integral 
relationship 

A 3 (e; E) = I dx (E 3 — 2xE 2 + x 3 ) £ (e, x) (13) 

In the periodic orbit theory, the correlation function (fTUj) can be expressed as a sum over 
classical periodic orbits [3J. The important energy scale in the system is that associated 
with the period of the shortest periodic orbit 
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For instance, in classically chaotic systems E max ~ A and for classically integrable systems 
-Emax ~ VeK For energies E E max , the levels are uncorrelated and one finds 

K(e 1 ,e 2 )-6(e 2 -E 1 ) (14) 
A 3 (e;E)~E/15 (15) 
S (e;E)~E (16) 

In the opposite limit, E 3> E max , the properties of spectral correlations are very different for 
the classically chaotic and classically integrable systems [lj. For the former, they are well 
known and are described by random matrix theory [7| and supersymmetric non-linear sigma 
model 8j. For the latter, it was believed that they lead to the saturation rigidity given by 

3 

A T(e;E) = J- 1 J2§ (17) 

3 i 

where Aj and 3} are the amplitudes and the periods of the periodic orbits and 2N is the 
dimension of phase space. 
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It turns out however, that the more precise formulae, up to the leading terms E mauX /E, 
are as follows: 

2 f(e 1 -e 2 )T j 



K c 



{El, £2) 



h N+l 



cos 



h 



A3 00 (s; E) 



A3 00 (e; 



h N - l A\ 



E 



COS 



£2 
ft 



(18) 
(19) 



sm 



ETj\ _-oo 

"2/T 1 : ^ 



.1- 



1 v°° t 2 



cos 



EI} 



where 



A3 00 (e; 25) 



h N 



9 4 2 



J J 



e-E) = 2Al 



(20) 



(21) 



In the above equations, the superscript "00" refers to "saturation behavior" and the overbar 
to averaging over the oscillations. Note that both Aj and Tj depend explicitly on the position 
of the center of the interval e ^> E. For instance, in a rectangular, with the aspect ratio of 



its sides L 2 /Li 



a, 



1/2 

asp ; 



one finds [3j that the periods are integers (representing the number 



of retracings) of irreducible cycles M = {Mi, M 2 } 

Tm. 



2h 



71 



eA 



2.,V2 , M 2 a -V2 
asp r J"2 ( -^ as P 



(22) 



where Mi and M 2 are coprime "winding numbers" of classical periodic orbits such that 

M 1 T 1 = M 2 T 2 (23) 
T12 being the periods of of motion along the sides L\ t 2. Expressions for Aj, and resulting 



formulae for the quantities of interest, can be found in Refs. [3J| and [lj. 

The key consequences of the above results are as follows. First, the "amplitude" of oscilla- 
tions around A^ (e; E) decays with the width of the interval E. Conversely, the "amplitude" 
of oscillations around S°° (e; E) does not decay with the increase of E; furthermore, this am- 
plitude is of the order of S°° (e; E). Second, the amplitudes of oscillations decreases rapidly 
with the period of periodic orbits. In a rectangle, for instance, A| oc T^ 1 and, using eq. 
( 1221) . it is easy to see that just a few terms with smallest winding numbers should dominate 
the sums in the above equations; this was indeed confirmed numerically [l|. 

To further appreciate these consequences, consider the contribution from a single har- 
monic only and compare the result with the known behavior of S°° (e; E) in completely 



uncorrelated system, in an almost rigid spectrum (Gaussian ensemble) and completely rigid 
spectrum (harmonic oscillator). It is convenient to consider the derivative <9£°° (e;E) /dE, 
for which we find |7[,[1[ 

Uncorrelated Nearly Rigid Rigid Oscillatory 

(integrable short range) (Gaussian ensembles) (harmonic oscillator) (integrable long range) 

dE K L (i ^ j -.short T~ S1U [ — 

where the summation is limited to the few shortest periodic orbits (and the corresponding 
harmonics). Clearly, depending on the interval width E, the oscillatory S°° (e;E) exhibits 
the range of behaviors, from uncorrelated to rigid. Moreover, <9£°° (e; E) / dE can become 
negative, implying a seemingly paradoxical result where the fluctuation of the number of 
levels decreases as the average interval width (and the mean number of levels) increases. 

Finally, because the frequencies of harmonics are incommensurate, S°° (e; E) ordinarily 
does not reach zero; it oscillates between E^ax ( £ > E) and S^ in (e; E), each typically of order 
of S°° (e; E) parametrically (for a particle in the box see |l(). However, if one can find a 
system where the shortest periodic orbit T dominates the sum, E°° (e; E) can be reduced, 
as per eqs. ()20| ]2l|) . to 

E«* (e; E) ~ 2S°° ( e; E) sin 2 

and can become effectively zero for E = 2hmiTQ 1 . The same effect would be also achieved 
if the periods of other periodic orbits are integer multiples of T . We found such a system 
in a modified Coulomb problem that we proceed to discuss below. 

III. MODIFIED COULOMB MODEL 

We consider a particle in the central potential 

V(r) = -- + 4 (24) 



Classically, the trajectory of the motion is given by 

V 



1 + ecos7 (9 — 9, 



19] 

(25) 
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where 



P= 2 -{e+€) (26) 



a \ 2m J 



Using the canonical action variables 







I r = -^L 2 + 2m(3 + aj^rry h = L (29) 



we can express the energy as 



ma 2 



2 (j r + + 2mf3 
and rewrite expressions for p and e as 



ma \I r + I l 

respectively. 

The frequencies of radial and angular motion are given by {13]: 



where the notation 



(30) 



P=— (31) 



^ = l^ = V^^ = 2 V2&i (32) 

de u r 

= ttt = — 33 
ale 7 



a 2 



was introduced. For any energy e, the motion is conditionally periodic except for the fol- 
lowing two circumstances. 

First, for the values of the angular momentum L such that 

M r 

7 = —— - rational (35) 
the motion becomes periodic with the periods of radial and angular motions related by 

T e = jT r or M e T e = M r T r = T M (36) 
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Here Tm is the period of an irreducible cycle M = {M r , Mg} (M r and Mg are coprime), 
whose orbital and angular winding numbers are, respectively, M r and Mg. In other words, 
the orbit closes for the first time after M r periods of radial motion and Mg periods of angular 
motion. 

Second, from 

e = <J=> I r = (37) 
the motion becomes circular for L such that 



e — 






e 





= ^2^^-^ (38) 
in which case M = {0, 1}. From eqs. (1281) and (I29p the corresponding frequency is found as 

4" r> = ^ (39) 

^ = vS (40) 

where u r is still given by (1321) but, since the distance from the center remains fixed, does 
not have the meaning of a radial frequency. 

It is very important to notice that u) r depends only on the energy e, and does not depend 
on the angular momentum L. As was already mentioned, at any energy e the conditionally 
periodic motion becomes periodic for such values of L that 7 is rational; these values, 
however, do not depend on e, except for the constraint 




7 = Wfr- -1=5 h 




2 

(cir) 



- 1 > 7 W (41) 



J \\\J |e| V2^p 

that follows from eqs. ( 1281) and (1401) . Consequently, the following picture of the periodic 
orbits emerges. In addition to circular orbits, whose period 

rp(dr) _ ^(cir) 

is given by eqs. (I3"2"|) and (HOT) and changes continuously as a function of energy, there are 
irreducible orbits such that 

M r = [Mei {cir) ] + i (43) 

where [] is the floor function and % are integers such that M r and Mg are coprime. These 
correspond to rational 7's (1351) and their period is given by 

T M = — M r = T r M r (44) 

UJr 



In view of inequality (|4TI) . new rational values of 7 become possible at discrete (quantized) 
values of e. In particular, the shortest periodic orbits 

M e = 1, 7 = M r = M r min + i (45) 

M min = jyrfr)] + ^ { < ^rfr) (45) 

are especially important. The key observation here is that as e increases (|e| decreases), the 

smaller values of M r become possible, with quantum jumps occurring for energies such that 
7 (cir) is 

integer. This fact will prove to be crucial in evaluation of the spectral rigidity and 
level number variance. Finally, for either type of the periodic orbit, it can be subsequently 
retraced with the period of nT-^, where n is the number of retracings. 



IV. SEMICLASSICAL SPECTRUM 



The quantum spectrum is given by 



2ma 2 ma 2 
£ P ,i = 2 5 (47) 

h 2 ^2p + 1 + y/(2l + l) 2 + 2m(3/h^j 2 ^ (p + V P + 2m P/ h2 ) 



which clearly follows from fl30|) via Born-Zommerfeld quantization of the action variables 

I r = h(p+l) ,I e = h(l + l) (48a) 



2 J \ 2 / 

and the second of eqs. (|47|) is the limit of large quantum numbers, p,l 3> 1, that in the 
semiclassical approximation used here. Further, we make yet another approximation that 
concerns with the fact that the standard Kepler problem (Bohr atom in the quantum limit) 
is a so called "supersymmetric" or "resonant" problem This is because the frequency of 
radial and angular motions coincide in the Kepler (so that the motion is periodic), which 
is indicative of extra symmetry in the problem, as well as additional conserved quantities - 
Runge-Lenz vector in the present case. In the Bohr atom, the latter is associated with the 
n 2 -fold degeneracy of the n's energy eigenstate. Conversely, in a "generic" (non-resonant) 
integrable system, the motion is ordinarily conditionally periodic, except for specific values 
of certain parameters upon which the motion may become periodic. Thus, in order to achieve 
the greatest possible difference with the standard Kepler problem, a large parameter (3 must 
be considered in f f24"|) . Accordingly, we assume that the condition /3/a 2 B ^> a /as holds, 
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where ae = h 2 /ma is the Bohr radius, that is m/3/h 2 ^> 1. This condition translates, as 
follows form (T4Tj) . to that of the quantum numbers p, I being limited from above which, 
combined with the semiclassical approximation, yields 



2m3 rad , , 

(49) 



Using (jl9"|) to expand eq. (JI7j) . we obtain 

- 2 - 2 2p^/2^3/¥ + l 
4/3 ' 4/3' 2m/3//l 2 ~ ^ ' ^ 2(3 

and, with a substitution, 



EpJ ra - ° + » jgV + ' 5 + ^ (50) 



(51) 

we find 

e Pii = 2p^2/3 + / 2 < 2/? (52) 
Classically, the condition corresponding to (J3SJ) would be 



h (= L),I r <£ (53) 

leading to 

2I r J2^ + I?; e 

£ ^ E «T 5 - E+ ^ (54) 

where 

e = 2/ rV ^+/ 2 <2/3 (55) 

and I r fi —>■ I r ,e/h are now dimensionless. 

Obviously, in both quantum and classical circumstances, e ~ — e in the zeroth order. The 
latter leads to simplified formulas for the frequencies since in this approximation 

\/2mp 7 

and, from eqs. (1401) and (J3 



7>7 (cir) = y^>l (57) 

as follows from (1531) . 

Due to the linear relationship between e and e in (15T?]) . spectral properties of the two are 
identical. Consequently, in what follows, it is spectrum (1321) that is studied numerically. 
Similarly to the rectangle, where ensemble averaging was understood in terms of variations 
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of the rectangle's aspect ratio jj, here ensemble averaging is understood in terms of the 
variations of (3. Flattening of the spectrum is achieved via ([2]) with 

£ 3/2 



<Af(e)> 



(58) 



which follows immediately from eqs. (CQ) and ( |52l) . This is equivalent to starting with the 
scaled Hamiltonian 

3/2 



(2I r VW + Ie) 



3^ 



(59) 



for which 



(Af{e sc )) = e sc (cN-sc) 
In what follows, we will drop subscript "sc." The frequencies are now given by 

1/3 



UJr 



where 




2(3 



(2(3 



> 7 (CJ ° = I ^ I » 1 



1/3 



(3e v ^75) 2/3 -2/ r v / 27J V3e 
obtained by solving ( |59j) for 1$. (For convenience, we carried over the notation 7^"")) 



(60) 



(61) 



V. LEVEL CORRELATION FUNCTION, SPECTRAL RIGIDITY AND LEVEL 
NUMBER VARIANCE 



We now turn to evaluation of the level correlation function ([18]) . Using the result obtained 
in Appendix, we find 

'2tt (ei-e 2 ) M r 



M 



U) r 

3eM r 



cos 



UJr 



(62) 



To reduce this to the sum on M r only, we notice that from (IBTj) and (TJ6 

M min = j 7 («r)j + j 



(63) 



Second, for each M r , there are [M r /^ c ' ir ^ possible values of M$, and, consequently, we 
obtain 

" M r 1 (2ir(s 1 -e 2 )M r 



M 



U) r 

3eM r 



1 



(cir) 



COS 



UJ r 



(64) 
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For very small (e± — 62), is possible to reduce this sum to an integral. Neglecting the 
difference between the function and its floor, including time-reverse of each orbit, and using 
eqs. ( IBTJj) and (IBTj) . we find 



K c 



E 

M r =M» n 
1 LU 2 r 

7T 3e7( cir ) 
5 (ei - £ 2 



2u; r 



3 e7 (cir) 



cos 



27T (5! - e 2 ) M r 



OJ,. 



where 



(65) 

(66) 
(67) 

(68) 

and T min is the period of the shortest periodic orbit. This is in complete analogy to the 
approximate form of the correlation function found for a rectangular box[l| where the 8- 
function term corresponds to the absence of level correlations and the second term the onset 
of thereof. 

Similarly, saturation spectral rigidity (|21[) is given by 



^ TV»1 



-1 



dx cos ((ei — e 2 ) x) 

lin 

sin ((ex -e 2 )/£) 
7t (ei - e 2 ) 



2P 



71^ 



E 



M r =M™ h 



7 



(cir) 



1 



7T- 



E 

M r =2 



7 



(cir) 



1 



(69) 



where E is now understood as interval width in the spectrum e. The second equality follows 
from the fact that the floor function in the sum automatically takes care of summation 
starting with M r min . 

Together, eqs. (f6"§|) . (I6"T|) and (|6"3"1) transparently predict quantum jumps in saturation 
level rigidity. As the energy increases, 7 C j r decreases and, as it takes on smaller integer 
values, a transition M™ 111 — > M™ 111 — 1 takes place leading to the jump in saturation rigidity. 
We observe such jumps in numerical simulations, discussed in the next Section. 

Finally, the level number variance is given by 



4^275 



7T- 



E 



M r =2 



M r 



7 



(cir) 



1 



sin 



fllEMr 



(70) 



assuring that S°° (e; E) = 2A^° (e; E). 
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FIG. 1: (e;E) vs. e for [3 = 3 x 10 6 . The thicker line is the numerical simulation while the 
thinner line is the analytical result given by eq. (|69p . with modification (|7ip . 

VI. NUMERICAL SIMULATIONS. 



We conduct numerical simulations on the spectrum (1521) for central values of (3 = 5 x 10 5 
and (3 = 1, 2, 3, 4, 5 x 10 6 . As previously mentioned, "ensemble averaging" is accomplished for 
each (3 by taking ~ 100 values of f3 around the central value. These /3's must be sufficiently 
close to the central value, as to eliminate the systematic dependence on f3, yet sufficiently far 
to ensure proper sampling. As a first step, we verified the Poisson (exponential) distribution 



for the nearest leve 



spacings, which should be the case for an integrable system without 



extra degeneracies |2|. 

The numerical result for vis-a-vis fl69|) are shown in Fig. 1 While fl69l) provides a 
wonderful fit to numerical data on the top plateau, it predicts much bigger jumps between 
the plateaus than is seen numerically. As a results, to fit lower plateaus we empirically 
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Z m (e,E) 




— I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - 1 ■— E 

1000 2000 3000 4000 5000 

FIG. 2: S°° (e;E) vs. E for = 3 x 10 6 and e = 5 x 10 5 . The thicker line is the numerical 
simulation while the thinner line is the analytical result given by eq. (|70p . with modification (|7ip 
(for this value of e, the correction is actually zero, as explained in text). 

added an appropriate empirical constant via substitution 

[M r /7< C "'>]-1 

+ E (71) 

n=l Z 

At this time we do not have a good explanation for this discrepancy. Further, the theoretical 
formula predicts additional noticeable jump further up in the spectrum, when ^ cir ^ = 3/2, 
which is not seen on the experimental curve. 

With the same modification ( 1711) in eq. ( 1701) . we plot E°° (e; E) as a function of E for two 
different values of e Clearly, while the periodicity is superbly predicted by the analytical 
expression, there is a striking difference in the intermediate structure in comparison with 
the numerical results. 

VII. CONCLUSIONS 

We predict analytically, and observe numerically, quantum Hall like jumps of the averaged 
saturation level rigidity in a modified Coulomb problem. These are explained semiclassically 
in terms of jumps in the winding numbers of the shortest periodic orbits as the position of 



M r 




M r 


o,(cir) 
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z 0> (e,E) 




, , , , 1 — i i i 1 i i i i 1 i 1 i i 1 i i L, 1 i E 

1000 2000 3000 4000 5000 

FIG. 3: Same as Fig. © for e = 2 x 10 5 . 

the interval center moves through the energy spectrum. Also, analytically and numerically, 
we predict sinusoidal oscillations of the saturation level number variance with the interval 
width. This is a striking result which indicates that while distribution of the levels on a 
interval varies greatly from sample to sample, the total number of levels in the interval 
may be nearly identical for certain values of the interval width. This is because all the 
higher harmonics are integer fractions of the shortest periodic orbit and add up coherently. 
The latter explains the difference with other systems, such as rectangular box, where the 
oscillations may be large but the variance doesn't reach a near-zero value. 
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APPENDIX A: AMPLITUDES OF PERIODIC ORBITS 

The amplitude of an irreducible cycle M = {Mi, M 2 } is given by[l 

2 27T 

m = T^u- dI M /dT M det {dui/dl k } m \ ( 1} 

uj (I M ) T m = 27rM (A2) 

This equation can be simplified by differentiating ( 1A2I) on T, which gives 

m So;, 9/,- „ , . „. 

3=1 

whereof 



" ' 7 " Tdet{^/a4} (A4) 

so thatQ 

M T M I- [w? (aj 2 /9/ 2 ) + {du x /dh)\ + ^u; 2 (cM/df, + du 2 /dh)\ M 1 ; 

which shows a universal dependence of v4 M on Tm- 

Here, I = {I r , Ig} and cj= {u r , ojq] = Vie (I), where e is given by (|59|) . Consequently, 
from ( 1A5I) . we find 

M ~ 3eT M 3eM r 

where ui r is given by fl60|) . 
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